IFUP-TH 44/92 



(N 



o 



> 

(N 



a 



Topological susceptibility and string tension in 
CP^-^ models S 



_ Massimo Campostrini, Paolo Rossi, and Ettore Vicari 

^^ I Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 

QQ ' Dipartimento di Fisica dell'Universita, F56126 Pisa, Italy. 

October 6th, 1992 



Abstract 



N-l 



o, 

(«_^ , We determined the topological susceptibility and the string tension in 

the lattice CP^~^ models for a wide range of values of A'', in particular for 
^^ ' A^ = 4, 10, 21,41. Quantitative agreement with the large- A'' predictions is 

5]^ ! found for the CP^° and the CP''° models. 

-)— > . 

^ , 1 INTRODUCTION 

(— I \ The most attractive feature of two dimensional CP™ ^^ models is their similarity 

with the Yang-Mills theories in four dimensions. Most properties of CP ~^ 
models have been obtained in the context of the 1/A^ expansion around the large- 
j^ ' A'^ saddle point solution. An alternative and more general non-perturbative 

JjHj \ approach is the simulation of the theory on the lattice. It is then important 

to perform Monte Carlo simulations of the lattice CP ~^ models to check the 
validity and the range of applicability of the 1/A^ expansion. 

The large- A^ expansion predicts an exponential area law behavior for suf- 
ficiently large Wilson loops, which implies confinement, with a string tension 
0{l/N), and absence of screening due to the dynamical matter fields. CP ~^ 
models have a non-trivial topological structure; at large A^ the topological sus- 
ceptibility turns out to be 0(1/A^). In this talk the features of CP ^^ models 
concerning confinement and topology are investigated. In order to study the ap- 
proach to the large-A^ asymptotic regime, we performed numerical simulations 
for a wide range of values of A^, in particular A^ = 4, 10, 21, 41. A more detailed 
discussion of these simulations and of their results can be found in Rcfs. p, 0, H| . 
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2 LATTICE FORMULATION 

We regularize the theory on the lattice by considering the following action: 

where Zn is an A'^-component complex scalar field, satisfying z„z„ = 1, and A„^^ 
is a U(l) gauge field. We also considered its tree Symanzik improved counterpart 
S*!^™ in order to test universality. 

The standard correlation length ^w is extracted from the long-distance be- 
havior of the zero space momentum correlation function Gp of two projector 
operators Pij{x) — Zi{x)zj{x). We define an alternative correlation length ^g 
from the second moment of the correlation function Gp. For N = 2 £,g/£,w — 1 
within 1% 0, while the large- A^ expansion predicts [|| 

Unlike ^w, which is a non- analytic function of 1/A^ around iV = oo, ^g can be 
expand in powers of 1/iV y|. 

3 TOPOLOGY 

The large- TV predictions concerning the topological susceptibility xt are [pi 

.2 1 f, 0.3801\ ^ / 1 \ ,„, 

and§ 

xta = 7^7 + (-ij,) . (4) 



Attn \N^/^^ 

Eq. ^ and Eq. (Q) are not in contradiction with each other due to Eq. (||) , but 
the first one should be testable at lower values of N according to the powers of 
N in the neglected terms. 

A troublesome point in the lattice simulation technique is the study of the 
topological properties and the determination of xt ■ While for large N the geo- 
metrical definition Xt lH is expected to reproduce the physical topological sus- 
ceptibility, at low N Xt could receive unphysical contributions from exceptional 
configurations, called dislocations. The dislocation contributions may either 
survive in the continuum limit, as it happens for the CP or 0(3) a model, or 
push the scaling region for Xt fo very large /3 values. 

Another approach, the field theoretical method, relies on a definition of 
topological charge density by a local polynomial in the lattice variables ql ■ The 



correlation at zero momentum of two q^ operators Xt i^ related to xt through 
the following equation: 

Xt^(/3) = a'ZiPfx, + a" A{(3){S{x)) 

+ P{m)+0{a^). (5) 

S{x) is the trace of the energy-momentum tensor, / is the identity operator. 
Z{P), P{/3), and A{P) are ultraviolet effects, since they originate from the ul- 
traviolet cutoff-dependent modes. The field theoretical method consists in mea- 
suring Xt^(/^)j evaluating Z{/3), A{f3) and P{P), and using Eq. (||) to extract 
Xt- 

A third method, called cooling method, measures Xt on an ensemble of 
configurations cooled by locally minimizing the action p2| . 

3.1 The CP^ model 

We determined Xt by performing simulations with both actions Sg and S'?''™ 
and by using different methods. In Fig. 1 the dimcnsionlcss quantity Xt^c '^^ 
plotted versus (,g- Data obtained by the geometrical method, indicated by the 
filled symbols in Fig. 1, violate universality, showing that A^ = 4 is not large 
enough to suppress the unphysical configurations contributing to xfi at least 
for ^ < 30. 

To apply the field theoretical method, we considered the following lattice 
operator: 

q'-i^) = -^ E ^M^Tr [P{x)A^P{x)A,P{x)] , (6) 

where A^ is a symmetrized version of the finite derivative. We neglected the 
contribution of the mixing with S{x): this assumption is supported by pertur- 
bative arguments. We obtained non-perturbative estimates of Z{/3 and P{/3) by 
using the "heating method" described in Ref. [Q. This method relies on the 
distinction between the fluctuations at ^ ~ a, contributing to the renormaliza- 
tions, and those at Z '^ ^ determining the relevant topological properties. Due to 
the critical slow down phenomenon, fluctuations at Z ~ a are soon thermalized, 
while, for large ^, the topological charge thermalization is much slower, allowing 
a direct determination of Z(f3) (when heating an instanton configuration) and 
P{/3) (when heating the flat configuration). 

An independent measure of xt is obtained using the cooling method. The 
results of both field theoretical and cooling methods are shown in Fig. 1. Scaling, 
universality, and good agreement between the two methods are observed. 

In conclusion, we quote for the CP'^ model Xt£,G — ^-^^ with an uncertainty 
of about 10%. 



3.2 The CP^ model 

We measured Xt by using the geometrical definition. Data were taken for ^q up 
to about 10 lattice units. Data for Sg show a slow approach to scaling, while for 
Sg^™ a better behavior is observed. For ^g the leading scaling violation term 
must be 0(ln ^/^^) when ^ ^ cxd Q . For the tree Symanzik improved actions the 
leading logarithmic corrections are absent, and scale violations are 0(^^^) |10| . 
In Fig. 2 we plot Xt Cg versus In ^g/Cg- Assuming that the scaling violation term 
proportional to In ^/^^ is already dominant in our range of correlation lengths, 
we extrapolated data of xf for the action Sg. We found XtCc — 0.0174(12), 
which is in agreement with the value of Xt^G obtained with the action S^^'"^. We 
then conclude that for the CP^ model Xt i^ ^ good estimator of the topological 
susceptibility. 

3.3 The CP2° and CP^° models 

As shown in Fig. 3 for both the CP^° and CP " models, data of Xt, obtained by 
using the geometrical method, are consistent with the large- A^ prediction (p|), 
whose results are indicated by the dashed lines. Fitting the data for the CP ° 
model, we found XtCc = 0.0076(3), to be compared with the value XtCc = 
0.00744 coming from Eq. M). However, the above result still disagrees with Eq. 
(||), which would require ^/Cw — 2/3, while we found ^q/Cw — 0-91 ^^ our 
simulations. 

In Fig. 4 we summarize our results by plotting Xt^a versus 1/-/V, showing 
how the topological susceptibility approaches the large- A^ asymptotic behavior 
(H), represented by the sohd line. 

4 CONFINEMENT 

The large- A^ prediction for the string tension cr is 

The string tension can be easily extracted by measuring the Creutz ratios 
defined by 

^(^'^)-l^M/(i?,T)W^(i?-l,T-l)- ^^> 

where W{R,T) are the Wilson loops constructed with the A^ field. In a 2-d 
finite lattice with periodic boundary conditions, the large abelian Wilson loops 
of a confining theory are subject to large finite size effects. For sufficiently large 
R the behavior of the Creutz ratios 'ri{R) = x{R,R), i-e. of those with equal 



arguments, should be |jl|: 



Tj{R) ^ a 



^ ,^-- 



(9) 



where L is the lattice size. To compare data from different lattices it is conve- 
nient to define a rescaled Creutz ratio 



7?,.(i?) = 7?(i?) 



L 



(10) 



In Fig. 5 we plot ti^{R)£,q versus the physical distance r — R/£,g for the CP^, 
CP^*', and CP***^ models. Starting from r ~ 2 the rescaled Creutz ratios show a 
clear plateau which is the evidence of the string tension. A similar behavior is 
observed for iV = 4. We do not see evidence of screening effects (at least up to 
r ~ 3^G, data for larger r were too noisy) confirming the picture coming from 
the 1/iV expansion. Data for the CP^" and CP ° show a good agreement with 
the large- iV prediction (0), represented by the dashed lines in Fig. 5. 

5 CONCLUSION 

With regard to topology and confinement, numerical results show a qualitative 
agreement with the continuum 1/A^ expansion for all values of N, while quanti- 
tative agreement is found for iV = 21 and for higher values of N. On the other 
hand, the approach to the large-iV asymptotic regime of quantities involving 
the mass gap appears very slow and the CP " model should be still outside the 
region where the complete mass spectrum predicted by the 1/N expansion could 
be observed M. The agreement in this sector of the theory is expected to be 
reached at very large N, because of the large coefficient in the effective expan- 
sion parameter Gn/N that can be extracted from a non-relativistic Schrodinger 
equation analysis of the linear confining potential [pj . 
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Figure captions 

Fig. 1: Xt£,G versus S,g for the CP^ model. We plot data of the geometrical 
method (for Sg "D" and ^l^" "O"), the field theoretical method (5g "O", S^^"^ 
"x"), and the cooling method (5g "D", S^^"" "o")- 

Fig. 2: Xt^G versus In^c/^g for the CP^ model. 

Fig. 3: xt^G for tfic CP^" and the CP'^" models. 

Fig. 4: XtS,G versus 1/A^. 

Fig. 5: r]r{R)^G versus r = R/^g, for N ^ 10, N ^ 21, and N = 41. For each 
N we show data at two values of (3. 



